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Stopping and regenerating a pair of single-photon pulses at adjacent locations in coherently pre¬ 
pared Rydberg atomic ensembles are significantly affected by their effective interaction mediated 
by Rydberg excitations, and the similar processes can differ notably from the one exhibiting the 
common Rydberg blockade as with the stationary propagation of multi-photon light beams in the 
same medium. Based on the complete dynamics, we reveal the detailed features in such processes 
by finding how the profiles of the involved quantum fields evolve in various situations. The findings 
help to determine the proper regimes for implementing photon-photon gates and photon transistors. 

In addition, we discuss the non-adiabatic corrections associated with quickly changing control fields, 
and illustrate a method that restores the photon pulses’ original amplitude during their retrieval 
unless they are heavily damped before storage. 


Since its first experimental observations [1|, 12| , the phe¬ 
nomenon of electromagnetically induced transparency 
(EIT) in ensembles of cold Rydberg atoms has attracted 
extensive studies. Different from the ordinary EIT Q, 
there exists van der Waals (vdW) or dipole-dipole inter¬ 
action between atoms, which modifies the dispersive and 
dissipative properties of light, as well as the photon cor¬ 
relations, in Rydberg medium [iSl. These interactions 
also indicate the possible formation of bound states 161 
and lattice structures [l7| of photons. More recent exper¬ 


iments 18-^ have demonstrated numerous interesting 
features related to Rydberg-atom-mediated interactions. 

The above mentioned long-range interactions make 
Rydberg EIT medium a promising candidate for quan¬ 
tum information processing with photons. A fascinating 
application of this kind is implementing a determinis¬ 
tic two-qubit gate, the building block of photonic quan¬ 
tum computing. Apart from the scenarios working with 
slowly traveling photons 23-221) a gate operation could 
be performed with stopped photons that are converted to 
Rydberg atomic excitations, so that the nonlocal interac¬ 
tion between the long-living Rydberg spinwaves achieves 
much more significant nonlinear phases than those gener- 
ated by the collisional interaction between atoms 28 , 23| ■ 
However, because the vdW interaction also b ring s about 
unwanted effects, the previous schemes [13, 311 in this 


category require more sophisticated technical steps, while 
it was not clear whether the simplest way to store and 
regenerate two photons (without extra procedures) prop¬ 
erly under vdW potential could realize a good gate oper¬ 
ation. On the other hand, other devices such as switches 


and transistors [33|, l3J should perform a different 


function that scatters more photons with a stored one. 
Stopping light in Rydberg medium and/or re-converting 
Rydberg excitations into light have been demonstrated 
experimentally [s^ - ls^ l. and the vdW interaction poten¬ 


tial in the used multi-photon beams is generally high and 
stable, leading to the well-known blockade effect. Yet the 
details of the corresponding processes for individual pho¬ 
ton pulses under their varying interaction, which are im¬ 
portant to the operations of the devices such as photon- 
photon gates and photon transistors, had not been fully 
understood. The previously unknown features in these 
processes will be illustrated in this work. 

In addition to the application concerns, the main pur¬ 
pose of our study is to provide a more realistic picture 
for the non-steady motions of interacting wave packets in 
Rydberg medium. To reflect the spatial distributions of 
photons and Rydberg excitations, we treat them as de¬ 
localized fields for their interactions and find their con¬ 
tinuous evolutions in space and time. This differs from 
most other related studies, which rely on a discontinuous 
blockade radius from point-point interaction potential. 

A typical example of our concerned problem is the 
setup of two pencil-shaped ensembles in Fig. 1. Here 
two single-photon pulses either travel along the same di¬ 
rection (co-propagation) or respectively enter the oppo¬ 
site tips of the ensembles (counter-propagation). The 
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FIG. 1: (a) Atomic level scheme. Here Ap = uieg — Wp, 
Ac = ojre — Wc (difference between the level gaps and the 
field frequencies), (b) Control field Rabi frequency Q.c{t) = 
ta,nh{tc — t)/Tc. Three values, Tc,i (solid), Tc ,2 (dash- 
dotted) and rc ,3 (dahed), with Tc,i < Tc ,2 < Tc, 3 , are given 
as the examples, (c) Geometry of the pulse propagations. 
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control fields with the Rabi frequency are being 

turned off with time, so that the pulses will be stopped 
and converted to Rydberg excitations that are close to 
each other. After the two Rydberg spinwave packets in¬ 
teract for a certain time, they will be turned on again 
to retrieve the photon pulses in the phase matched di¬ 
rections. This process can be modeled by the dynam¬ 
ical evolutions of the electromagnetic fields f;(x, t) of 
the two photons (^ = 1 and 2) with the kinetic Hamil¬ 
tonian Hp/h = —icJ2iJ {x)dz£i{x.), the induced 
polarization fields Pi (x, t) = '/Nag^ (x, t) and Rydberg 
spinwave fields Si (x, t) = yfN (x, t ), where Upu = 
A high atomic ensemble density N enhances 
the photon-atom coupling (with the constant g) seen in 
the Hamiltonian Haf/h = — f dx{g\/]V£^(x)Pi(x) + 
^c(i)Sl(x)Pi(x) + II.c.} + f dxApP/Pi(x), which de¬ 
scribes the atomic level scheme. A narrow band photon 
propagates with negligible absorption under the EIT con¬ 
dition Ap -I- Ac = 0 [Ap/c is defined in Fig. 1(a)]. How¬ 
ever, under the spinwave interaction with the Hamilto¬ 
nian = J dxj dx'Sl(x)S 2 (x')A(x — x')S 2 (x')Si(x), 
where A(x — x') = Cq/\x — x'\^ is the vdW poten¬ 
tial, the condition will be violated by a shift V of the 
levels |r) in the relevant atoms. The consequent dis¬ 
sipation from populating the levels |e) that decay at 
the rate 7 can be depicted by a stochastic Hamiltonian 
Hdis/h = dx{Cj{x,t)Pi{x) - H.c.} involving 

the quantum noise fields 0 (x, t) of the reservoirs [i^ . 

It is sufficient to find the mean values Oi{x,t) = 
{Oi{x,t)) of the quantum fields Oi{x,t), and their evo¬ 
lutions under the total Hamiltonian H = Hp -|- Haf -f 
Hint + Hdis are determined by the following equations 
(see Appendix A and we set = 1 here): 

dt£i[x,t) + cdz£i{x,t) = igVNPi{x,t)] ( 1 ) 

dtPi{x,t) = + iAp)Pi{x,t)+iVll{t)Si{x,t) 

+ ig'/N£i{x,t)-, (2) 

dtSi{x,t) = -iVi{x,t)Si{x,t) +incit)Pi(x,t), (3) 

where 

= y dx'|—^^|S' 3 _i(x',t)p (4) 

for ^ = 1 and 2. Using Eqs. (ED and one can expand 
the right-hand side of Eq. © as 

,— . . q^N d . . 

zgV^P,(x,t) = -^^-£,(x,t) 

g‘^N£i{x,t) d 1 gN dr 1 ,d 

Hc(<) ilcit) dt^n*(t)^dt~^ ^ 


Then, in our concerned regime realizing slow light 
( 5 ^ A /1 Hep 1 ), the time derivative on the left-hand 
side of Eq. o will be absorbed into the leading term 
in the above equation, reducing Eq. (1) to the one only 
with a spatial derivative. The discretization of this rear¬ 
ranged form of Eq. (1) is used to get the spatial profile of 
£i at a specified moment, given the values of Pi induced 
by the field £i at the preceding moment of iteration with 
the discretized Eqs.© and ©. We apply a second or¬ 
der Runge-Kutta method in the iterative steps toward 
the quantities (f;. Pi, Si and Vi) all over the space-time 
grids. For simplicity, only symmetric propagations in two 
ensembles will be discussed below, but the generalization 
to different pulse motions is straightforward. 

As a time-dependent boundary condition for pulse 
evolution, the electromagnetic field of the input pho¬ 
ton pulses is supposed to have the profile £lp{p,t) = 
Jo( 2 r'oip/d) at the entries of the ensem¬ 
bles, where is the maximum of the photons’ Rabi 
frequency Sip = g£, and tp and Tp are the time scales in¬ 
dicating the peak arrival and pulse duration, respectively, 
while Jo{x) is the Bessel function of order zero with its 
first zero point t'oi. The effects of a gradually increasing 
interaction can be best seen from stopping two counter- 
propagating pulses. In reality, storage of photons may 
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FIG. 2: Dynamics of the photon pulses (a, c) and the induced 
spinwave in the unit (b, d), which counter-propagate 

in two 100-/im long ensembles of ^’^Rb atoms subject to a 
control field of Tc = 1 gs, tc = 100 gs, and = 1.5 x 27r 
MHz in Fig.l(b). The field profiles are on the ensemble axis, 
and obtained from the numerics with the iteration step size 
0.002 gs in the time direction and 0.02 gm in the spatial 
direction. Here, \g) = 5Si/2, |e) = 5Pi/2, and |r) = 1005'i/2, 
with Ce ~ —2.3 x 10® GHz gnP and 7 = 5.75 x 27r MHz. We 
set Ap = —Ac = 57 in (a, b) and Ap = —Ac = —57 in (c, 
d). The photons are Gaussian pulses with = 0.01 MHz, 
tp = 10.0 gs and Tp = 5.0 gs. The system parameters are 
A = 2 X 10^® cm“®, a = 6 gm and d = 2 gm. The insets are 
the corresponding imaginary (solid) and real (dashed) parts 
of the normalized susceptibility to the integers. 
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FIG. 3: (al) Realized spinwave profiles (on the ensemble axis) from stopping two counter-propagating pulses in two 300-/rm long 
ensembles, given the different ensemble separations and a switching of Tc = 10 /rs, = 45 /rs. (bl) Corresponding potentials 
V(z) the standstill spinwave packets in (al) create on the other ensemble. (a2) and (b2) show the similar quantities for the 
co-propagating setups, while the control field is turned off at tc = 30 /rs. Ensemble length for co-propagation is flexible as 
long as the together pulses can be contained inside. A nearly homogeneous potential around 0.02 MHz is realized with the 
largest a in (b2). (cl) Spinwave profiles from counter-propagating pulses stopped at different switching speeds (a = 10 /rm); 
(c2) similar profiles for a co-propagating setup, (dl) Example of counter-propagating spinwave dynamics, given a = 10 /im, 
Tc = 10 ^s, and tc = 40 /rs. (d2) Example of co-propagating spinwave dynamics, which differs from (dl) by tc = 24 fis. These 
examples demonstrate that different propagation geometry can achieve similar results, but the co-propagagation setups enjoy 
the advantage of much shorter ensembles. The common parameters are = 2 x 27r MHz, — 0.01 MHz, Tp = 7 /rs, tp = 12 
US, N = 2 X 10*^® cm“®, and d = 2 /rm. The iteration step sizes in the numerical simulatious are the same as those in Fig. 2. 


fail as shown in Fig. 2, which illustrates the associated 
dynamical evolutions for the photons with the opposite- 
sign detuning values. There a positively detuned photon 
will go through the medium, though it is damped. If the 
photon’s detuning changes the sign, it will not survive the 
travel. Such difference can be explained with the suscep¬ 
tibility x(wp), defined as P{tOp) = \/]Vx(wp)fIp(a;p), for 
a single-frequency wave under a constant interaction po¬ 
tential R; see the insets of Fig. 2. In the former situation 
the attractive potential V “pulls” the central frequency 
component initially under the FIT condition away from 
the absorption peak, but it “pushes” the corresponding 
component in the latter toward the peak. Then the pos¬ 
itively detuned photons will virtually see two-level sys¬ 
tems when they get closer, but the same medium becomes 
opaque to the negatively detuned ones. These effects are 
more significant for a higher detuning Ap, as the gap be¬ 
tween the absorption peak and the two-level regime nar¬ 
rows down accordingly. The scenario of negative detun¬ 
ings is most suitable for implementing a photon transis¬ 
tor. In the presence of many photons at the same time, a 
significant interaction potential can immediately “push” 
the system across the absorption peak to the two-level 
regime on the other side, fitting the Rydberg blockade 
model (super-atom) irrespective of the photons’ detun¬ 
ing signs which well describes the stationary propaga¬ 
tion of light in Rydberg medium. By contrast with the 
phenomenological models, our first-principle calculations 
provide the understandings of what happen to individual 
photons under their constantly varying interaction and a 
time-dependent flc(t) simultaneously. 

In order to realize their storage, the photons, as well 


as the control field, should be better to be resonantly 
coupled to the energy levels (Ap = Ac = 0) so that the 
farther positions of absorption peaks leave more room for 
the frequency components not to enter the above men¬ 
tioned regimes. Certainly one can take a stronger control 
field to widen the FIT window, but the practice sup¬ 
presses the induced spinwave, as seen from the relation 
S ~ g'/N £ jflc for a single-frequency wave without in¬ 
teraction. Accordingly the photon pulses should have 
a sufficiently long duration Tp (corresponding to a nar¬ 
rower bandwidth 6ujp = l/rp) to fit into the limited FIT 
window. Such spatially extending photon pulses will be 
mapped to Rydberg spinwaves distributing over the en¬ 
sembles, after the control field is turned off. When the 
process happens to two pulses counter-propagating re¬ 
spectively in close ensembles, the mutual interaction can 
diminish their fronts before they are stopped side by side 
(the interaction between their back portions is still in¬ 
significant at this time), leaving a spinwave packet of 
asymmetric shape as in Fig. 3(al). The potential V{z) 
acting on other ensemble can therefore be extremely in¬ 
homogeneous as shown in Fig. 3(bl). For an ideal pho¬ 
tonic gate operation, each point on the stored wave pack¬ 
ets should be under the same potential value V, so that 
the wavefunction will gain a uniform phase. Thus one 
has to face a trade-off between the magnitude and the 
uniformity of the achieved potential V{z). 

Unlike most other previous researches, we adopt the 
potential field in Fq. o, which is generated by the real¬ 
time spinwave field S'(x, t), for the the mutual interaction 
between Rydberg excitations. The evolution an input 
pulse undergoes under this potential is very different from 













4 


that under the interaction with a single-point Rydberg 
atom. The interaction potential value V (x, t) is simply 
the degree of deviating the EIT condition Ap -|- = 

0 at any space-time point (x, t), when the control field 
is on. The consequent space-time dependent dissipation 
can be captured by our simulations directly based on 
Eqs.((T])-(l3l). More examples of such pulse dissipation are 
given in Appendix B. The calculations also specify the 
proper regime of larger ensemble separations, in which 
an approximately homogeneous potential can be realized 
for gate operation; see Eig. 3(b2). In this regime beyond 
certain ensemble separation, co-propagation and counter¬ 
propagation actually make little difference. 

The pulse motion in ensemble is largely influenced by 
the control field. Given the Rabi frequencies in Fig. 1(b), 
the parameters tc and Tc determine where a pulse will be 
stopped. A fast switching with small Tc can immediately 
stop the pulses at any moment. One question is whether 
the non-adiabatic corrections [s^, connected with terms 
containing the time derivatives of on the right-hand 
side of Eq. could lead to considerable pulse losses. 
This can be clarified by our numerical calculations with 
arbitrary which go beyond the previously studied 

adiabatic passages [38l-l40j. The results in Figs. 3(cl) and 
3(c2) indicate unremarkable difference between the dis¬ 
parate switching speeds. Further calculations manifest 
that, in the ensembles of high densities, the non-adiabatic 
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FIG. 4: Dynamical evolutions of photon pulses throughout 
the whole “write-in” and “read-out” process. The storage 
period is the same as in Fig. 3(dl), with = 2 x 27r MHz, 

Tc ~ 10 fis and a = 10 /rm. The control fields for retrieval 
have their Rabi frequencies Q,c'^^{t) as the time-reversals of 
the functions in Fig. 1(b). We take = 1 x 27r MHz 

in (a) and = 3 x 27r MHz in (b), while the switch- 

on occurs at Tc = 0.1 /rs in both cases, (c) Profiles of the 
retrieved photons for different ensemble separations, where 
— 5 X 277 MHz and Tc = 0.1 jj.s. The dashed curve is 
that of the input photons for comparison, (d) Profiles of the 
photons regenerated with different switching speeds as shown 
at the exits, given = 2 x 27r MHz and a = 10 /rm. 


effects are overshadowed by those from the pulse interac¬ 
tion; see Appendix C. 

Under any circumstance the photon absorption is in¬ 
evitable due to pulse interaction, limited EIT window, 
non-adiabatic effects and others. We here demonstrate a 
mechanism that compensates for the previous dissipa¬ 
tions through photon regeneration. In Fig. 4(a) the 
control field in the “read-out” is weaker than the con¬ 
trol field for the “write-in”. The amplitude of the re¬ 
generated photons is less than that of the input pho¬ 
tons, and they exit the medium at a lower group ve¬ 
locity (flatter contour) as well. The retrieval control 
field in Fig. 4(b) is, on the other hand, stronger than 
the one used for stopping the photons, so that the 
re-converted photon pulses become brighter and leave 
the ensembles at a higher group velocity (steeper con¬ 
tour). Near the entries and exits of this counter¬ 
propagation setup (almost no interaction exists there), 
the pulses can be regarded as dark-state polariton field 
'I'(x, t) = cosi?(t)f (x, t) — sini?(t)S'(x, t), where d(t) = 
arctan gy/N/^c(,t). The electromagnetic field of 

the pulses therefore takes the form cos'd(tm)4>(x, U„) 
at the beginning and cos'd (tout )'f'(^Uout) in the end. 
To maintain their amplitudes, the dissipation that re¬ 
duces the field 'F(x, tm) to 'if {x, tout) can be partly off¬ 
set by the factor cosd(t), though the retrieved pulses 
are deformed (such deformation does not affect carrying 
quantum information). The best enhancement render¬ 
ing cost)(tout) = 1 for the output photons demands an 
infinite U°"‘, implying the impossibility of restoring the 
amplitudes of the severely damped pulses with a realis¬ 
tic U°”‘(t) [See Fig. 4(c)]. Because pulses’ dissipation is 
primarily decided by the duration and intensity of their 
interaction, a quick retrieval outdoes a slow one, con¬ 
trary to what happens to the storage period; compare 
Fig. 4(d) with Fig. 3(cl). A slow “write-in” and a fast 
“read-out” are thus favored to our concerned processes. 

In summary, we have presented the first study of the 
details in stopping and regenerating individual photon 
pulses under their mutual interaction in Rydberg EIT 
medium. Much more complex and varied than those of 
the stationary propagation of multi-photon beams, sev¬ 
eral important features in the processes originate from 
the ever-changing interactions which lead to space-time 
dependent dissipation. Our simulations provide more re¬ 
alistic pictures of the continuous photon pulse evolutions 
in space and time, dispensing with a boundary of dis¬ 
continuity as the blockade radius in phenomenological 
models. Apart from clarifying the regimes to implement 
photon transistors and gates, and ascertaining the in¬ 
significant role of the non-adiabatic corrections in high 
density ensembles, we show that the amplitudes of the 
weakened photon pulses can be enhanced by a stronger 
retrieval control field. These understandings could be 
valuable guide for the relevant experimental researches. 
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A. Dynamical equations of quantum fields 

In our concerned problem, the photons with the electromagnetic fields fi(x, t) are co-propagating or counter- 
propagating in two ensembles, having their kinetic Hamiltonian (here we give the form of co-propagation) 

Hp/h = -ic'^ f dx£j{x)dz£i{x). (A-1) 

1=1 

Together with the control field of the Rabi frequency they will induce the polarization fields Pi{x,t) and the 

Rydberg spinwave fields Si (x, t) in the medium, as described by the many-body version of the atomic-level-scheme 
Hamiltonian 


Haf/h = 


Z=1 


x)P/(x) -\-nc{t)Sj{x.)Pi{:K.) -\-H.c.} 


2 

E 

* 


dxApP/(x)P/(x), 


(A-2) 


where g is atom-field coupling constant, N the atomic density of the ensembles, and Ap the detuning of the photons. 
The Hamiltonian for the interaction between the Rydberg spinwave fields takes the form 


Hint = da dx'^J(x)^|(x')A(x - x')52(x')5i(x). 


(A-3) 


Meanwhile, as the result of the interaction, the population of the intermediate levels spontaneously decaying at the 
rate 7 gives rise to the photon dissipation, and we describe it by the stochastic Hamiltonian 

Hdis/fi = f dx{Cj{x,t)Pi{x) - H.c.} (A-4) 

1=1 

in terms of the coupling between the polarization fields and quantum noise fields, where the random-variable noise 
operators satisfy [^i(x, f), ^)„(x', t')] = 6imS{x — x')6{t — t'). An advantage of this stochastic Hamiltonian approach 
is that the physical process can be modeled by the formal unitary evolution operator U{t, 0) = Texp{—i P(T)dT} 
with H(t) = Haf + Hint + Hp + Hdis , which naturally leads to the Heisenberg-Langevin equations of the relevant 
quantum fields as shown below. 

In deriving the dynamical equations, we write dP;(x, t) = Q(x,t)dt, which satisfy the Ito’s rules (a generalization 
from those in [sH'): 

dBi{x,t)dBi{x,t) = 0, dB\{x,t)dB\{x,t) = 0, 

dB\ (x, t)dBi{x, t) = 0, dBi{x, t)dBj (x, t) = dt. (A--5) 

Then an increment of the polarization fields, for example, will be found as 
dPi{x, t) = U^{t + dt, t)Pi{x., t)U(t -\- dt, t) — Pi{x, t) 

= -^[Pi{x,t), [Haf + Hint P Hp)dt-ih^/^ J dx'[dB{x',t)P^{x',t) - dB\x',t)Pi{x',t))] 

-I- 7 y (ix'(2P'^(x',t)P(x,<)P(x',t) — P(x,t)P'l'(x',t)P(x',f) — P'^(x',t)P(x',f)P(x,t)), (A-- 6 ) 

where the above Ito’s rules have been applied to the expansion up to the second order of —iHdis{t)dt. In this way 
the exact dynamical equations for the three types of quantum fields are obtained as follows {h= 1 ): 

dt£i{x,t) P cdz£i{x,t) = igy/NPi{x,t)-, (A-7) 


dtPi{x,t) = -{-iPiS.p)Pi{yi,t)piVtl{t)Si{x,t)Pig'/N£i{x,t)-^/^li{x,t)-, 


(A-8) 


7 




FIG. S-1: Comparison between the potential from a point particle (a) and that due to a stopped distributing spinwave (b) 
for different ensemble separations a. Here the potentials originate from a vdW potential ~ in Eq. dm. The profiles 

in (b) from stopping two counter-propagating pulses are achieved with a control field of = 2.0 x 27r MHz, Tc = 8/is, and 

tc = 40 fis. When these obtained spinwaves shrink to one point on the ensemble axis, they give the potentials in (a). The 
photons resonantly coupled to the energy levels have = O.OlMHz, tp = 12.0 fis and duration Tp = 7.0 fis. The atomic level 
scheme is the same as the one in Figs. 2-4 of the main text, and the atomic density of the 0.335-mm long ensembles with their 
diameters d = 2 pm is N = 2 x 10^® cm~®. 


dtSi{x.,t) = -iVi{x,t)Si{:x:,t) + (A-9) 

where 

Vi(x,i)= y dx'A(x - x')^]_,(x',t)53_i(x',t) (A-10) 

for Z = 1 and 2. As in other photon storage scenarios [s3| . we will replace the field operators in the above dynamical 
equations by their mean values to see how their profiles evolve under given conditions. The dynamical equations for 
the field profiles in Eqs. (l)-(3) of the main text are thus obtained by averaging out the quantum noise operators, 
(6(x,t)) = 0. 


B. More examples of pulse evolution under realistic interaction 

The interaction potential a pulse experiences in the concerned processes is determined by the spatial extension of 
the other Rydberg spinwave packet, as indicated by Eq. (lA-lOl) or Eq. (4) of the main text, which comes from the full 
dynamics discussed in the above. This leads to a significant difference from the interaction with a point Rydberg atom. 
We illustrate the fact with the comparison in Fig. S-1. The finally obtained realistic potential profiles in the illustrated 
situations correspond to those of the standstill Rydberg spinwaves. Though the extending spinwave packets lower the 
potential magnitudes, the potential values distribute less abruptly over the ensembles. Such potentials due to wave 
packets instead of point particles exist for any type of interaction. For instance, the dipole-dipole interaction potentials 
between the stored wave packets should also be evaluated in this way, and are found to be more inhomogeneous than 



FIG. S-2: Realized dipole-dipole potentials from A(x — x') ~ Gs/lx — x'p in Eq. (IA-3II . These potentials are created by 
applying a static electric field parallel to the ensemble axis, after the photons are stopped and converted to Rydberg spinwave 
packets. The atomic level scheme is the same as the one in Figs. 2-4 of the main text, to have ICsj « 6.65 x 10® MHz pm®. 
The other parameters of the 0.3-mm long ensembles are the same as those in Fig. S-1. Due to the large ensemble separations, 
the accompanying vdW potentials almost take no effect, resulting in the same spinwave profiles as seen in (a). 












those of vdW interaction; see Fig . S-2. The dipole-dipole interaction is applied in another photonic gate proposal 
based on stored photon pulses S3| . More realistic considerations should be therefore included in quantum information 
processing by the similar methods. 

One of the meaningful consequences of such realistic interaction potentials is that they lead to different dissipation 
rates over the pulses. In Ref. [s^ about another topic related to Rydberg FIT, the interaction between pulses is 
also modeled as the one between extending wave packets, but their dissipation is assumed to be with a global rate 
independent of the pulse distributions. The exact details of pulse dissipation should be clarified by approaching the 
processes with the complete dynamics of the quantum field profiles. Here we present two more examples to show how 
such space-time dependent dissipation manifest in photon storage. The first one in Fig. S-3 is about the even closer 
ensembles in Fig. 3(al) of the main text. The interaction between the pulses in these setups becomes stronger. Since 
the counter-propagating pulses have spatial extensions, the interaction between their fronts is intense before they are 
stopped in the medium. An extreme result is that half of a pulse will be absorbed. Negligible decay of the Rydberg 
levels after the control field is turned off maintains stable asymmetric spinwave profiles during the storage period. 
The other example in Fig. S-4 illustrates the dynamical evolutions of two co-propagating photon pulses and their 
associated spinwaves before they are stopped together. Here we use a narrower pulse bandwidth (corresponding to 
longer pulse durations Tp) so that the pulse dissipation is almost completely due to their interaction. This example 
clearly demonstrates the significant effects of pulse interaction. Due to a longer interaction time, only a small portion 
of the initially induced spinwaves will remain in the end. 
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FIG. S-3: (a) Created spinwave profiles for various relatively small ensemble separations, (b) Dynamical evolution of the 

associated spinwave field on the axis of ensembles. Here the ensemble separation is a = 6 pm. The control field is switch off at 
a speed of Tc = 10 ps, and ensemble length is 0.3 mm. The other parameters are the same as Fig. S-1. 
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FIG. S-4: (a) Dynamics of two co-propagating photon pulses during their storage, (b) Dynamics of the associated spinwaves. 

Here we take a = 10 pm, tp = 30.0 ps, Tp = 18.0 ps, Tc = lOps, and tc = 80 ps. The other parameters are the same as those 
as in Fig. S-1. 


C. Further discussion on non-adiabatic corrections 

We now take a further look at the effects from a fast-changing Rabi frequency Hc(t) of control field. To neglect the 
photon losses from these effects, the systems should evolve along adiabatic passages during storage and retrieval of 
photons. The condition for realizing such adiabatic passages can be formulated in terms of a change rate of the mixing 
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FIG. S-5: (a) Relations between the storage efficiency rj = \S^(ts)/{tin)\ and the control field switch speed indicated 

by Tc, where and \S^{tin)\ are, respectively, the peak value of the stored and initial spinwave. These results for 

different atomic densities N are about the storage processes without pulse interaction, (b) Corresponding relations for two 
co-propagating pulses in the ensembles separated by a distance a = 13 /rm. In both situations we take Tp = 10 fis, tp = 20 fis, 
and tc = 40 /rs. The ensemble length is 6 mm. The other parameters are the same as those in Fig. S-1. 



angle = arctang'/N/fldt) as t? <C \\E± — Eq\, where \E± — Eq\ is the energy difference between bright-state 
and dark-state polaritons (see, e.g. [s^). The change rate, which takes the form 



gVN 

g^N + ndt) 


ndt), 


(C-l) 


scales inversely with the square root of the atomic density N in our concerned regime realizing slow light 
{gVN /|flc(t)| !)■ In ensembles of high density N, insignificant non-adiabatic effects are therefore expected. This is 

seen from Fig. S-5(a) about the processes without the interaction between pulses, in which the storage efficiency goes 
up with the density and the difference between fast and slow switching of a control field is not so large. Meanwhile, in 
the presence of the pulse interaction as shown in Fig. S-5(b), the tendency of storage efficiency with atomic density 
will be totally reversed. The interaction between spinwaves becomes stronger in an ensemble of higher density, because 
the spinwave profiles S'(x, f) = \/]V(d'gr-)(x, f) are magnified by the square root of N. Accordingly more significant 
dissipation will exist in the ensemble. These results demonstrate that the dissipation from non-adiabatic corrections 
is hardly comparable with that due to pulse interaction. Stopping and regenerating photons in the setups can be done 
with quickly switching control field too. 
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